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Use of a cDNA microarray to 
analyse gene expression 
patterns in human cancer 

Joseph DeRisi'% Lolita Penland : & 
Patrick O. Brown 2 (Group ] ); 
Michael L. Burner 3 *, Paul S. Meltzer 3 , 
Michael Ray 3 , Yidong Chen 3 , Yan A. Su 3 & 
Jeffrey M. Trent 3 (Group 2) 



The development and progression of cancer 1 " 3 and 
the experimental reversal of tumorigenicity 4 * 5 are 
accompanied by complex changes in patterns of 
gene expression. Microarrays of cDNA provide a 
powerful tool for studying these complex phenom- 
ena 6 " 8 . The tumorigenic properties of a human 
melanoma cell line, UACC-903, can be suppressed 
by introduction of a normal human chromosome 6 t 
resulting in a reduction of growth rate, restoration 
of contact inhibition, and suppression of both soft 
agar clonogenicity and tumorigenicity in nude 
mice 4 - 5 - 9 . We used a high density microarray of 
1,161 DNA elements to search for differences in 
gene expression associated with tumour suppres- 
sion in this system. Fluorescent probes for 
hybridization were derived from two sources of cel- 
lular mRNA [UACC-903 and UACC-903(+6)] which 
were labelled with different fluors to provide a direct 
and internally controlled comparison of the mRNA 
levels corresponding to each arrayed gene. The flu- 
orescence signals representing hybridization to 
each arrayed gene were analysed to determine the 
relative abundance in the two samples of mRNAs 
corresponding to each gene. Previously unrecog- 
nized alterations in the expression of specific genes 
provide leads for further investigation of the genet- 
ic basis of the tumorigenic phenotype of these cells. 

DNA microarrays, containing I J6I total elements, 
including 870 different cDNAs and controls'** 1 1 {see 
Methods), were printed rohotically onto a glass micro- 
scope slide in four quadrants covering an area of about 
1 enr (Fig. ] ). We prepared fluorescent cDNA probes 
using total poly (A) + mRNA from UACC-903 cells and 
UACC-903(+6) cells by labelling with a green and red 
fluor, respectively. A mixture of the two flourescently 
labelled probes was hybridized to the DNA microarray. 
This comparative hybridization method, coupled with 
the doping of synthetic standards and an estimation of 
statistically significant deviation for local background 
variance allowed a direct and quantitative comparison 
of the relative abundance of individual DNA sequences 
in this complex sample'* - *. We added a set of synthetic 
poly (A) + -tailed 'mRNAs* to the purified mRNA from 
each cell line as internal standards to assist in quantita- 
tion and estimation of experimental variation intro- 
duced during labelling and reading. Targets 
complementary to these standards were included, in 
duplicate, on the microarray. Based on these standards, 
mRNA species comprising 1:10,000 of the mass of the 
poly ( A) + RNA could readily be detected. 

In a representative two-colour fluorescent scan of all 



spond to genes preferentially expressed in the tumori- 
genic UACC-903 cell line, and the reddish spots corre- 
spond to genes preferentially expressed in the 
non-tumorigenic UACC-903(+6) cell line. Genes 
expressed at approximately equal levels in the two cell 
lines appear yellow or brown. A portion of the array at 
higher magnification highlights the diverse pattern of 
differential expression observed (Fig. 2b). In Fig. 2c rec- 
tangles corresponding to specific array elements are 
coloured to reproduce the hue and intensity of the fluo- 
rescent signal at each element. The hybridization signals 
from a duplicated set of genes are shown juxtaposed, to 
illustrate the reproducibility of the hybridization signals 
for each gene. 

To address the possibility that an apparent difference 
in expression might result from experimental variables 
unrelated to the difference in chromosomal composi- 
tion between the two cell lines, we examined the vari- 
ance in expression for 90 'housekeeping 1 genes. We 
selected these genes based on the assumption that they 
would not be differentially expressed between the two 
cell lines. The averaged red/green ratio for this subset of 
genes was 1.13. The averaged red/green ratio for the set 
of five internal standards was 0.97 ( m = 10). The vari- 
ability in the expression level of the housekeeping genes 
probably overestimates the experimental variability in 
measuring differential expression. As a conservative stan- 
dard, an absolute fluorescent signal (red or green) with 
an intensity greater than that observed at the contr I 
array elements containing total human genomic DNA 
was considered to represent specific hybridization. Gene- 
specific hybridization was therefore only considered sig- 
nificantly different between samples if the following two 
criteria were met: i) the signal intensity (green or red) 
exceeded this threshold; and ii) the logarithm of the 
red/green fluorescence signal ratio differed by >3 S.D. 
from the mean logarithm of this ratio for the 'house- 
keeping gene panel (that is, ratios <0.52 or >2.4). 

By these criteria, mRNA levels for 15/870 (1.7%) genes 
were significantly diminished, while the mRNA levels 
for 63/870 (7.3%) genes were significantly increased in 
association with suppression of tumorigenicity by intro- 
duction of chromosome 6. To test the reliability of 
microarray hybridization results in identifying differen- 
tially expressed genes, we analysed 16 genes by north- 
ern analysis, in each case, the results of northern analysis 
corroborated the differential gene expression identified 
by microarray hybridization (Fig. 3). 

Significant differences in expression between these 
two cell lines identified several genes as candidates for 
determining features of the tumorigenic phenotype of 
the melanoma cells. For example, among the genes 
detected with significantly higher expression (>10-fold) 
in the tumorigenic cells was the human brown locus pro- 
tein (TRP1 /melanoma antigen gp75). This is the most 
abundant glycoprotein in mcianocytic cells and a critical 
melanosome membrane protein 1213 . Additionally, its 
expression is reduced when melanoma cell lines are 
induced to differentiate by treatment with HMBA I2J3 . 
Also expressed at a significantly higher level was a spliced 
variant of the mRNA encoding myelin IM.P/DM20. This 
is widely expressed in neural crest derived cells in early 
development and has been suggested to play a role in 
cell-cell signaling during development M . 

.'. i. c :..'.j.ir/. ; .! -.. ...^.\'Alev- 
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^^^^^^^^^^^^^^^™T™~"~T" . . nNA —Med onto a ooiv-tvsine coaled slide. The DNA is stained with a DNA-specific fluorescent 
Bg. 1 Properties of cDNA microarrays. a. A "j-orescen ^^^^^S^mil up to 10.000 spots/2.54 X 7.62 cm microscope stide. b. Erfi- 
dye. YOYO. The cenler-to-center ,dT)- to arrays in the absence of competitor produces strong 

cient blocking o« ^'"^ oCS such as the EST TM827 shown. Rhodamine-labelled cDNA (red) Irom the UACC-903 cell 

hybridization to immobilaed poly (dAC as well as t0 ^^^ u ^^^'^t ai H = , ota , human) . similarly, hybridcation with fluorescein-labelled 
Je hybridized in the presence of poly «*)• Mocker sh «m human SSI and atsomecDNA elements (presumed to con- 
Coll DMA in the absence of ^^ & ^ , ^^eWM ™ the UACC-903 cell line produces little it any signal a. these 
tain highly repeated sequences, such as R23416). DNA Tne abse nce of signal at some cDNA locations following 

locafons when hybndized in the presence^ of exces i l ^^Zi^J^^S^^ cDNA targets do no. contribute s,gni.ican. hybridiza- 
UACC-903 cDNA hybridizat.ons also indlM,e V™,£!j^^^ was carried out with 4 print heads, spaced to fit into 4 adja- 

£nS,^^^ 

target types in each of the resulting quadrants is shown. 



els were elevated by the addition of a normal chromo- 
some 6 ( 1 7 genes) are known lo be activated by IFN-y. a 
cardinal proinflammatory cytokine that, among other 
activities, induces expression of the gene products ot the 
MHC class II locus. For example, the mRNA encoding 
monocyte chemotactic protein 1 ( MCAF/MCP1 K a 
cytokine that induces monocyte chemotaxis and activa- 
tion IM< \ was more than 10-fold less abundant in the 
tumorigenic cell line. In the skin, M PCI is critical in the 
regulation of cutaneous monocyte trafficking 1 *" 1 ', and 
elevated expression plays a role in suppression of tumour 
growth and metastasis 1 ^ 21 . The mechanism by which 
these interferon-y regulated genes are induced in UACC- 
903 cells by transfer of a normal chromosome 6 remains 
to be determined, it is worth noting, however, that the 
interferon-y receptor gene is localized to the distal long 
arm of human chromosome 6. 

Finallv, several genes that showed > 10-fold higher 
expression in the suppressed UACC-903(+6) cells have 
previously been recognized in other models of tumour 
suppression. Most notably, there was elevated expres- 
sion of the mRNA encoding WAFI (p21), a key media- 
tor of tumour suppression by p53 (ref. 18). The p21 
protein had previously been identified as a melan ma 
differentiation-associated antigen (termed mda-6) ,9 * :o . 
In melanoma cell lines suppressed for metastasis by the 
introduction of chromosome 6, expression ot WAFI 
(p21 ) mRNA and protein correlates inversely with 
metastatic potential 20 . 



These results provide a wide view of the diverse sys- 
tems that are altered in this model system of tumori- 
genicity, and focus attention on specific gene products 
and pathways that may be of particular importance in 
this tumour type. 

Our ability to classify human cancers in a way that 
reflects the underlying molecular pathology or that 
anticipates their potential for progression or response 
to treatment, remains primitive. Using cDNA microar- 
rays to define alterations in gene expression associated 
with a specific cancer may be an efficient way to uncov- 
er clues to the specific molecular derangements that con- 
tribute to its pathogenesis and thus identify potential 
targets for therapeutic intervention. Moreover, recogni- 
tion of pathognomonic alterations in gene expression 
might provide a basis for improved diagnosis and mol- 
ecular classification of cancers and thus allow selection of 
the most appropriate therapeutic strategies. 

Public databases of human expressed gene sequences 
contain partial sequences of at least 40,000 different 
human genes 11 , and efforts to develop a human tran- 
script map have developed rapidly 31 . Based on the high 
yield of information obtained using an array of < 1,000 
different genes, a more comprehensive survey of gene 
expression patterns, using a more complete array of 
human genes, will likely provide a rich source of new 
and useful insights into human biology and a deeper 
understanding^ the gene pathways involved in the 
pathogenesis of cancer and other diseases. 
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Fig. 2 DNA microarray analysis of changes in gene expression between the tumorigenic cell line. UACC-903. and its non-tumorigenic derivative UACC903< + 6> 
derived by introduction of a normal chromosome 6. a. A ratio image of the results of simultaneous hybridization of Rhodaminel 10-labelted cDNA ioreeni 
from UACC-903 and Cy3-»abel>ed cDNA (orange-red) from UACC-903<*6) to a microarray. To produce this image, the 

orescent probe were combined as the appropriate colour channels in a single image. Arrows indicate the location within the array oUhe corresoondino 
genes analysed by northern blotting (F,g. 3). b. A magnified image of the area of the array boxed in white in (a), c. Magnif.ed image of trJee cOrWSe^ 
0* %JZ W V n J a) ' re ^ em ' n9 the t C ? NAS f f0r ™ '^ AF/MCP ' 1 W ratio >10,; centre. fracUn (r/g ratio 1.04); andnght. ^^ZS^^W^ 
0.2) [see Fig. 3]. d. simplified representat.on of ratio hybndization results. Quantitative fluorescence intensity data is extracted from each array taraet The aver- 
age target colour ratio determines the hue of each box and the average intensity determines the brightness of each box. In this image, the ordeToYthe boxes 
corresponds to their or.g.nal order in the microliter plate from which they were printed. Duplicate printings of the same plate can be examined side bv side 
as in the first two rows shown here, to assess reproducibility of the hybridization results (see text). Numbered arrows indicate the location within the array cor^ 
responding to genes analysed by northern blotting in Fig. 3. ' 



Methods 

'Generation of microarrays, hybridization, scanning. The 
preparation ol coated microscope .slide?* and .subsequent robot- 
ic priming ol DNA was carried mil in a manner .similar to that 
described', briefly, pre-cleaned glass slides were treated with 
poly-l. -lysine solution (Sigma) to form an adhesive surface lor 
printing. PCM products,' purified by ethanol purification, were 
resuspended in 3.\ SSC. A custom built arraying robot picked 
up and deposited small volumes (-5 nanoliters) of DNA onto 
the slides. Alter printing, the slides were washed in a 0.2% SOS 
solution. The remaining hound DNA was denatured by sub- 
merging the slides in 95 "C distilled water for 2 min followed by 
a brief wash with ethanol. DNA was W crosslinked to the 
slides (Stratagem* St rata I inker, M) ml). 1o prevent non-specific 
probe binding, the slides were blocked by rinsing in a solution 
of 70 mM succinic anhydride dissolved in O.I M boric acid pH 
8.0, containing 35% I - met hyl-2-pyrrolidi none (Aldrich). 
Additional protocols and parts list pertaining to microarrav 
fabrication can he. obtained from http://cmgm.stanford.edu/ 
pbrown. 

Puritied. labelled cDNA was resuspended in 1 1 pi of 3.5xSSC 
containing 4 ug i>f poly (dA)* DNA, 2.5 p.g fi. ru/ilKNA. 4 Mgof 
human Cot I DNA (Ciibco IMI.). and 0.3 ul of It)'!,, SDS. Prior to 
hybridization, the solution was boiled for 2 min then allowed to 
cool to room temperature. Hybridization was carried out at 
f»2 °C lor -U h in a water hath. Prior to scanning, slides were 
washed in 2* SS< !. 0.2"... SI >S for 5 min and 0.:>: SS( ! for I min. 

local laser microscope built bv S. Smith with software written 



by N. Xiv. A separate scan, using the appropriate excitation line, 
was done for each of the two lluorophores used. Data was col- 
lected at a maximum resolution of y microns/pixel with 12 hits 
of depth. 

Probe preparation and labelling. UNA was extracted from cells 
using the Tria/ol reagent (IT1 Inc.l, following the manufactur- 
ers directions, cl )NA probes were synthesized from singly oligo 
dT-selcctcd (Pharmacia) mKNA pools. Huorescently labelled 
cDNA was prepared from mKNA by oligo dT-primed polymer- 
ization using Superscript II reverse transcriptase (LTI Inc.). 
The pool of nucleotides in the labelling reaction was 0.5 mM 
dGTP. dATP and dCIP and 0.2 mM dTTP. Fluorescent 
nucleotides. Rhodaminc lit) dUTP (Perkin Rimer Cetus) or 
Cy3 dUTP ( Amersham), were present at 0.1 mM. Probes were 
purified by get chromatography (liioSpin 6/HioKad) and 
ethanol precipitation. 

Selection of cDNA elements and generation of control tem- 
plates. Synthetic cDNAs were prepared by cloning random 
/inwHI and Hmdlll ended fragments of E, colt DNA in the vec- 
tor pSPo-l poly (A)' (Promega). linearizing isolated plasmid 
DNA with EuM\ and synthesizing poly (A)* tailed UNA com- 
plementary to the insert from the resident SPd promoter 
(Piomcga). Prior to use, the synthesized UNAs were selected on 
oligo dT cellulose. The largest group of cDN As consisted of 674 
cDNA clones from the iNlfi arraved normalized infant brain 
l.i;iar» . liie.se lio;;. , -,•.„.,, seated to *wry NIB 

library member that corresponded to a named gene according 
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MARCKS 




Collagenase 



MCAF/MCP-1 



Fig. 3 Nonnern hybridization suDstantiai- 
ing the consistency of the cDNA rrucroar- 
ray results. Corresponding locations within 
the cDNA microarray illustrated in Fig. 2a 
are provided for 1) Waf-Vp21; 2) MARCKS\ 
3) collagenase; 4) MCAF/MCP-1; 5) n-1- 
antichymotrypsin; and 6) /* -actin. The sig- 
nal detected by a radio-labelled 0-acfm 
probe represents a control for loading vari- 
ance, with a red/green ratio observed on 
the cDNA microarray (Fig. 2a,c) for fLactin 
of 1.04. 
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io the UniGcnc EST clustering sys- 
|cm :u: p hc j^-onj largest group of 
clones consisted of IH3 sequenced cDNA 
clones general ed by .subtraction of c DNA 
from l he chromosome-ft suppressed 
non-lumorigenic UACC-903 (+6) cell 
line with cDNA from its parental tumori- 
genic cell line UACC-903 (ref. 9). 
Approximately 100 additional genes 
(total 870 genes arrayed) were obtained 
from KST libraries on the basis of their 
a- 1 -Antichymotrypsin expression pa I tern 1 1 ismic specific, and so 
on). Kach array included the following 
hybridization controls: plasmid vector, 
lambda, OX 1 74 phage, total human DNA, 
. human Coll DNA, and poly (A}\ The 

P~ Acttn synthetic Mandards used for normaliza- 

tion of signals in each wavelength were 
also arrayed. Controls were included in 
each quadrant of the array to assess the reproducibility of the 
hvbridizaiion signal. Two plates of cONA clones (derived Irom 
the UACOM3 subtracted library) were also arrayed in dupli- 
cate. Pideliiy of the Unigene array relative io dbKS T was tested 
by sequencing of a random sample of II clones used lor 
microarrav construction. All sequences were identical with the 



corresponding dhKST entries. Additionally, each muuuiMWii 
cDNA from I he l'ACC-903 subtracted library was scquciucd. \ 
listing of cl>NAs comprising this microarray iWiuh weie 
derived from the Unigene and 'housekeeping' panel »..in he 
obtained from http://\vuw.nih.gov/l)IU/LCC/AKKAV/'e.\pn.- 
html. 

Northern blot analysis. Total UNA, 10 ug per lane, was elec- 
trophoresed in 1.2% agarose- formaldehyde gels and transterrcd 
onto nylon membrane (Hybond-N\ Amersham) by capillary 
blotting overnight. For DNA probes insert fragments from the 
Scares I Nib cDNA library 1 " were obtained by vector PCR for 
p2l, MAUCKS. a-l-aniichymoirypsin and (i- act in. Probes for 
fibroblast collagenase and MCAF/Ml'P- 1 were isolated from a 
UACC-V03(+M enriched cDNA library 4 with all probes 
labelled bv random priming. Fillers were washed to a strin- 
gency of 0. Ix SSC at 12 °C for 20 min. 

Web sites, http://cmgm.stanford.edu/pbiinvn for protocols and 
parts list pertaining to microarray fabrication. 
http://www.nchgr.nih.gtw/DIK/l-C(;/AKKAY/expn.html for a 
listing of cDNAs comprising this microarray which were 
derived from the Unigene and 'housekeeping' panel. 
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